What are the best routes for us to use for driving home tonight in rush hour traffic? 
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We show that the capacity of a complex network that models a city street grid to support congested 
traffic can be optimized by using routes that collectively minimize the maximum ratio of betweenness 
to capacity in any link. Networks with a heterogeneous distribution of link capacities and with a 
heterogeneous transport load are considered. We find that overall traffic congestion and average 
travel times can be significantly reduced by a judicious use of slower, smaller capacity links. 
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Choosing the best route for driving home from work is 
a task that many of us face everyday. Knowing the best 
choice is also important to urban planners who design 
city transportation networks. To answer this, and other 
important questions, a vast amount of research has been 
devoted in recent years to the analysis and optimization 
of transport on complex networks 111, l 2|, [ a , |4| . la, l a, M , 
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If we would collectively choose to use the best routes, 
then the traffic congestion and our commute times could 
be reduced. When traveling in a city, the best route 
to use is presumably the one that takes the least time 
to go from origin to destination. If there is little or no 
traffic, this route is simply the so-called shortest-path 
route. On the complex network that describes the city's 
streets, the shortest-path route is defined as the path for 
which the sum of the weights of each of the links along 
the path is minimal. The weight of a link can be defined 
as the time it takes to travel from one end of the link to 
the other when traffic is light. Generally, these weights 
are inversely related to the transport capacities of the 
links (streets that can handle a lot of traffic are traveled 
faster). On the other hand, when traffic is heavy some 
links may become congested, and then the quickest route 
to the destination may be a longer one involving smaller 
capacity links. It is not obvious, however, exactly at what 
level of congestion the use of alternative routes becomes 
advantageous, or how much the improvements in network 
transport capacity and average travel time amount to. 
In this paper we present methods to help answer these 
important questions and demonstrate their effectiveness 
within the framework of a simple model that captures the 
most important characteristics of urban street networks. 

Specifically, we use the complex network model intro- 
duced by Barthelemy and Flammini in [3]. This model 
has been shown to produce planar geometric networks 
with an exponential distribution of the node distances 
from the origin and other characteristics similar to the 
actual street networks of various cities. We use this 
model to demonstrate the effectiveness of our methods, 
rather than a particular realistic example, in order to 
make a statistical analysis by taking ensemble averages 



over many example city networks. Unlike previous stud- 
ies d lU [13, lH, 14, 15 1, which assumed that traffic is 
limited by node congestion, here we assume that the ca- 
pacity of the network is limited by the amount of traffic 
each link can support. This important and realistic dif- 
ference requires a nontrivial variation in our methods. As 
we have shown previously 12 , 3, l3| , when transport is 
limited by node congestion, the transport capacity of the 
network is maximized by using a set of routes that min- 
imize the maximum betweenness of the nodes. We also 
introduced an algorithm for finding such a set of routes, 
and demonstrated that, for a number of commonly stud- 
ied network topologies, it finds routes for which the ca- 
pacity, at least, scales optimally with system size. In this 
case, the transport capacity of the network is maximized 
by minimizing the maximum betweenness to capacity ra- 
tio of links. Here we also consider heterogeneous link 
capacities and uneven traffic demands between the vari- 
ous pairs of origin and destination nodes. In particular, 
we study the case of a rush hour traffic burst emanat- 
ing from a central location. To obtain our results, we 
use a variant of our previous routing optimization algo- 
rithm. Furthermore, we prove a formula that allows the 
quick computation of the average of the sum along the 
path of any link-related quantities and use this formula 
to compute average travel times. 

In our network transport dynamics, particles are as- 
sumed to travel along the network according to static 
routing rules. To demonstrate the importance of our re- 
sults, we consider two different types of routing. One is 
the "natural" shortest path (SP) routing computed with 
link weights in inverse proportion to their respective ca- 
pacities. The other one is the "optimal" routing (OR) for 
congested traffic resulting from the application of a rout- 
ing optimization algorithm which finds the set of routes 
that maximizes the transport capacity of the network. 
Note that the OR is the collective optimum that occurs 
if everyone uses the best routes for collective results. As 
such, it describes an important limit of collective behav- 
ior. This contrasts with the goal of the multitude of 
traffic studies that seek to optimize results individually 
through a learning process within the framework of evo- 
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lutionary game theory [J]. 

Here we assume that each outgoing hnk of a node has 
a separate "first in, first out" (FIFO) particle queue. 
Transport on the network proceeds in discrete time steps 
and is driven by inserting new particles at the nodes. The 
average number of new particles inserted per time step 
at node i with destination at node j is rij and we denote 



Each new particle is inserted at the end of 



the appropriate queue at its node of origin, namely the 
queue corresponding to the first link it has to traverse on 
its way to destination. The transport capacity of a link 
(i, j) is Cij, defined as the average number of particles 
transported per time step assuming an infinite number 
of particles in its queue. We use a stochastic sequential 
updating of the positions of particles that erases the cor- 
relations between particle arrival times and ensures that 
both the arrival and the delivery of particles are Poisson 
processes IJ] . The load of the network is defined as the 



average (r) over all nodes on the network of . The net- 
work transport capacity is the critical value (r)^ of the 
load above which the network becomes jammed. 

Routing on a network is characterized by the set of 
probabilities for a particle with destination t cur- 
rently at node i to be forwarded to node j. The between- 
ness of a node i with respect to a source node s and 
a destination node t is defined as the sum of the prob- 
abilities of all paths between s and t that pass through 
i. Node betweenness can be computed in terms of the 
probabilities pj;*'^ ^,^2^. The betweenness of link 
going from i to j is defined by [24 1 
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The average number of particles crossing link {i, j) in the 
course of a time step is given by 



N 
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To account for an uneven traffic demand pattern be- 
tween the nodes while still using a single parameter to 
characterize the load of the network it is convenient to 
write Tij = C (r) pij , where (r) is the average number of 
particles per node generated in the course of a time step 
and are nonnegative demand weights {pa 0) . If the 
weights are normalized such that Y^^j^i Pij = N{N — 1) 
we find 



{r) 
N - 1 
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and, using Eq. [2l the flow of particles through link [i,]) 
becomes 
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FIG. 1: (Color online) Average travel time (in time steps) 
vs. network load (in particles per time step per node) for a 
typical network realization with — 200 nodes. Results are 
shown for shortest-path routing, SP (blue dashed lines), and 
for optimal routing, OR (red solid lines), at three different 
high capacity link densities, (a) Ph — 0, (b) Ph = 0.3, and (c) 
Ph = 0.7. 



with the total weighted betweenness of link (i, j) given 
by 



N 
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Avoiding congestion means Wij < Cij for every link 
(i,j). Consequently, maximum transport capacity is 
achieved when the highest betweenness-to-capacity ra- 
tio on the network {B /C)max is minimized and is given 
by (r)^ = {N — \) / {B /C)max- The minimization of 
{B/C)max can be achieved by various methods includ- 
ing, but not limited to, simulated annealing and extremal 
optimization. 

The average travel time can be computed from link 
betweennesses assuming the arrival and the delivery of 
particles at every node are Poisson processes. We will 
now prove a general formula for calculating averages over 
the entire set of routes considered on a network and then 
apply this formula to calculate the travel time. Let Qij 
be any quantity associated with the links of a network. 
To calculate the average over all particle routes of the 
sum of Qij along the route it is convenient to write the 
betweenness in terms of the probabilities for complete 
routes. Let 7r„(s, t) be the ordered set of all nodes along 
the n-th distinct route between s and t (including s but 
excluding t) and Pn{s,t) be the probability for a particle 
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FIG. 2: (Color online) (a) Ensemble averages of the network 
average and network maximum betweenness-to-capacity ra- 
tios, {{B /C)avg) and {{B/C)max) respectively, vs. average 
number of nodes at = 0.3 for networks with an aver- 
age of 250 nodes. Results are for {{B/C)l^g) (black solid), 

{{B/C)Z1) (red dotted), ((B/C)^L) (green dashed), and 

(^(B/C)^ax') (blue dot-dashed), (b) Exponents of the power 
law fits for the same quantities as functions of P^. 



to be routed along that route. The number of distinct 
routes between s and t is Af{s,t). Then, 

A^(s,t) 

^ij'*^ = X! X! ^ik ^],next{k\n,s,t} (6) 

n=l kGiT„{s,t) 

where next{k\n, s,t) is the successor of node k in 
7r„(s, t) U {t}. Let us now compute the quantity {T,Q)i^vg^ 
defined by 
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By substituting Eq. [5] into Eq. [3 we find 

{^Q)avll' = E ^'"(*'*) E <3fc,nea;t(fcln,s,t)- (8) 
n=l fce7r„(s,t) 

The inner sum on the right-hand side of Eq.[S]is the sum 
of Qij along the n-th route and (S(3)it,g^ is its average 
over all routes between s and t. 

The network-wide average of {^Q^avg' weighted by the 
elements of the traffic demand matrix rgt is then 



avq 



(9) 



where Bij are defined in Eq. [5] and we used the fact that 
the sum of all Pst equals iV(iV — 1). The derivation above 
parallels the one presented in [1^, but the difference in 
the way one accounts for the links along a route in Eqs.[6] 
and [Has opposed to accounting for nodes is non-trivial. 

The average time, measured in steps, required for a 
particle to traverse link (i, j) in the absence of additional 
delays is 



1 + 
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where {qij ) is the average length of the queue associated 
with this link and dj is the link capacity. Since both 
the arrival and the delivery of particles at every queue 
are well approximated by Poisson processes, the average 
queue length is given by [25] 



Wi. 



(11) 



By using Eqs.fTOlandfTTjinlQjwe obtain the average travel 
time as a function of load 
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Additional time delays associated with traveling along 
the links may also be included in the calculation of Tavg 
by using Eq. [5] 

To find the set of routes that maximizes trans por t ca- 
pacity we use a variant of extremal optimization 2^, 27 1 . 
Specifically, we modified an earlier algorithm that has 
been shown to converge to near-optimal routing in the 
case of both geometric and small-world networks if 
traffic is limited by the processing capacities of the 
nodes ij, |lj, ll5| . The original algorithm aims at min- 
imizing the maximum betweenness of any node on the 
network. The variant employed here minimizes the max- 
imum betweenness-to-capacity ratio of any link on the 
network. It is, therefore, capable of maximizing the 
transport capacity of networks with traffic limited by ar- 
bitrary link capacities. The algorithm works as follows: 

1. Assign an initial weight equal to one to every link 
and compute the shortest paths between all pairs 
of nodes and the betweenness of every link. 

2. Find the link that has the highest betweenness-to- 
capacity ratio {B / C)max a-ud increase its weight by 
one. 

3. Recompute the shortest paths and the between- 
nesses. Go back to step 2. 

The computational complexity of the algorithm is 
0{N'^ log TV) (it is 0(iV2 log TV) for one iteration and 
requires 0{N) iterations to converge), which is the same 
as that of the variant that applies to transport limited by 
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node capacity [l4|. The algorithm, therefore, provides an 
approximate solution in only polynomial time to what is 
known to be an NP-hard problem As we will show, 
the approximate solutions we find, at least, scale opti- 
mally with network size A'^. Furthernore, because the 
resulting optimal routing is static, it only needs to be 
calculated once. Thus, the algorithm is useful even for 
very large networks. 

For simplicity, we consider networks with a binary dis- 
tribution of link capacities. Then, "streets" are low ca- 
pacity links and "highways" are high capacity links. Each 
network realization consists of a network topology and a 
set of capacities, generated by choosing any link to be 
a highway with probability Ph. Networks are generated 
using the algorithm introduced by Barthelemy and Flam- 
mini [T^ , which is based on the idea of links (highways 
or streets) growing gradually towards population centers. 
All links grow at the same constant rate towards the cen- 
troid of their adjacent population centers, and new pop- 
ulation centers are generated at constant average rate. 
When a link reaches the centroid of its adjacent popula- 
tion centers it splits into separate links directed towards 
each center. At the end, the nodes of the network will be 
the population centers and the points where links have 
been split. This model has been shown to generate net- 
works with characteristics similar to those of real-world 
street networks, particularly if the distances of the pop- 
ulation centers with respect to the city center are ex- 
ponentially distributed. These similar characteristics in- 
clude the distribution of the node degrees, as well as size 
and shape distributions of the areas delimited by streets. 
Using this model rather than a particular real-world ex- 
ample allows the calculation of ensemble averages over 
thousands of model city networks. An ensemble of net- 
work realizations is characterized by the number of cen- 
ters generated per time step Rc, the rate of growth of 
the links d, the total number of centers Nc, the ratio 
between the capacity of a highway and that of a street 
Ch/Cs, and the probability Ph for a link to be a high- 
way. Note that the number of nodes on these networks 
can only be controlled on average, through Nc and the 
ratio Rc /d. Without loss of generality, the average dis- 
tance of the population centers from the center of the city 
is assumed to be unity. All results presented below have 
been computed for Rc — 0.1, d — 0.001, and Ch/Cs — 8. 
Due to the high value of Ch/Cs, shortest-path routing 
will typically forgo the use of lower capacity links unless 
they are strictly needed to achieve connectivity. 

Figure [1] shows a comparison between the SP and OR 
average travel times for one network with uniform traffic 
demand at three different values of Ph- The total travel 
time along a link is taken to be the waiting time in its 
queue plus an additional delay equal to the inverse of the 
link's capacity. The waiting time in queues accounts for 
congestion, while the inverse of a link's capacity accounts 
for its travel time in the limit of low traffic. Note that 
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FIG. 3: (Color online) Ensemble averages of the network 
average and network maximum betweenness-to-capacity ra- 
tios, {{B /C)avg) and {{B/C)max) respectively, as functions 
of Ph for networks with an average of 250 nodes. Results 
are for {{B/C)^Q (black solid), {{B/C)°^g) (red dotted), 
((B/C)^^,) (green dashed), and ((B/C)^!) (blue dot- 
dashed) . 



routing is explicitly optimized to maximize the transport 
capacity and not to minimize the average travel time. 
The results shown are typical for networks with an av- 
erage of about 200 nodes. The load above which OR 
outperforms SP is in this case about 80% of the criti- 
cal load under SP routing. This may be regarded as the 
threshold for congestion when SP routing is used. Opti- 
mization increases the transport capacity of this network 
by a factor of about 1.7 if all routes have the same capac- 
ity, while for values of Ph between 0.2 and 0.8 the factor 
of improvement increases to more than 2. Note also the 
decrease of the average travel time at low loads with in- 
creasing Ph regardless of the type of routing. This is due 
to the increased average capacity of the links. 

In Fig. [5] results are presented for the ensemble av- 
erages of {B/C)avg and {B/C)max for networks with 
(N) « 250. All error bars represent 2cr estimates. The 
subscripts avg and max denote the average and respec- 
tively maximum over all links of a network with a given 
link capacity configuration. Averaging over an ensemble 
of network realizations is denoted by angular brackets. 
This is an average over an ensemble of network topolo- 
gies and over an ensemble of link capacity configurations 
for each network topology. An ensemble of network real- 
izations is characterized by the average number of nodes 
(N) and by the probability Ph for a link to be a highway. 

Figure mja) shows a log-log plot of {{B/C)avg) and 
{{B/C)max) for both SP and OR at Ph = 0.3. Since the 
lines are nearly straight, the ensemble averages of these 
quantities scale with average network size as a power 
law. Note that any decrease in {{B/C)max) can only 
be obtained at the expense of an increase in {{B /C)avg) 
since avoiding congestion along the shortest path means 
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taking longer routes that contribute to the betweenness 
of more links. With OR, the slopes of {{B/C)avg) and 
{{B /C)rnax) are essentially the same. Thus, the capac- 
ity of the routes we find, at least, scales optimally with 
system size. This behavior is similar to that of the node 
betweenness when transport is limited by node process- 
ing capacity. However, the finite size effects are stronger 
in the current case and the error bars under estimate 
the true error, since they are calculated assuming that 
the values of [B/C)avg and [B/C)max for the various 
network realizations are normally distributed while our 
simulations show they are not. Figure [2l^b) shows the ex- 
ponents of the power law scaling of the four quantities as 
functions of Pu- Note that the exponents for {{B /C)^^) 
and {[B/C)'^^) are essentially equal over the whole 
range of Ph which argues in favor of the optimality of 
routing. Note also that the exponent for {{B/C)^^^) 
exhibits a dip around P^, — 0.9. This is an interesting 
feature, indicating that the SP routing works unusually 
well when there is a small but nonzero concentration of 
low capacity links. 

Figure[3]shows < B/C > versus Ph for {N) « 250. The 
error bars in this figure also represent 2a estimates. Note 
that {{B/C)avg) varies monotonically between P^ = 
and P/i = 1 corresponding to all streets and all highways, 
respectively, while {{B / C)max) exhibits a midrange max- 
imum in the case of SP routing that corresponds to a dip 
of the network transport capacity. This type of behavior 
is due to the fact that SP routing forgoes the use of low 
capacity links as long as they are not strictly needed to 
achieve connectivity (just as we tend to do in real life), 
which increases congestion on highways. On the other 
hand, by optimally using all links, it is possible to avoid 
this phenomenon. 

A comparison between the SP and OR average travel 
times for one network with uneven traffic demand is pre- 
sented in Fig.m Specifically, we look at a "rush hour traf- 
fic burst" with particles originating from the innermost 
10% of the nodes that we use to model a "downtown". 
Their destinations are chosen uniformly from among all 
nodes with the only constraint being that the node of 
origin must be different from the destination. This is 
the same network that was used for Fig. [1] and again the 
results are representative for networks with an average 
of about 200 nodes. Note the lower maximum capacity 
of the network, which is due to the uneven distribution 
of traffic. Nevertheless, a judicious use of the low ca- 
pacity links again results in a significant increase of the 
transport capacity. The factor by which transport capac- 
ity is increased may be even higher than in the case of 
uniform traffic demand. This result is particularly impor- 
tant since traffic congestion usually develops is situations 
of uneven traffic demand. 

In summary, we have studied how to optimize collec- 
tive transport on planar networks that model a city street 
grid. Traffic on these networks is limited by link conges- 
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FIG. 4: (Color online) Average travel time (in time steps) 
vs. network load (in particles per time step per node) for a 
"rush hour traffic burst" with particles originating from the 
innermost 10% of the nodes on a typical network realization 
with N = 200 nodes. Results are shown for shortest-path 
routing, SP (blue dashed lines), and for optimal routing, OR 
(red solid lines), at three different high capacity link densities, 
(a) Ph = 0, (b) Ph = 0.3, and (c) Ph = 0.7. 



tion. The network model accounts for many topological 
features seen in real life urban street networks and al- 
lows for heterogeneous link capacities and uneven traffic 
demands. The question posed in the title is answered 
in statistical terms, as we show that a judicious use of 
all links, including low capacity ones, can significantly 
alleviate traffic congestion. The results we present show 
that the average factor by which transport capacity can 
be increased varies with network size as a power law, and 
that the routes we found are ones that allow the capacity 
to, at least, scale optimally with system size. Moreover, 
the average travel times on congested networks can be 
significantly decreased. We determine the network load 
at which using the low capacity links becomes advanta- 
geous from the point of view of the average travel time. 
Most importantly, we show that significant improvements 
in traffic routing can be achieved not only in the case 
of uniform traffic demand, but also in situations of ex- 
tremely uneven traffic. Our results are important not 
only for quantifying the routing improvements that can 
be achieved on existing networks, but also for the design 
of future street networks. For example, numerical exper- 
imentation with optimal traffic routing on a model of a 
real city street network can pinpoint the links for which 
an increase in capacity will result in the highest network 
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transport capacity. 
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